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1 Introduction 



Consider a ferromagnetic model, without quenched randomness, evolving from a disor- 
dered initial state, according to some dynamics at fixed temperature T. In the high- 
temperature paramagnetic phase (T > Tc), the system relaxes exponentially to equilib- 
rium. At equilibrium, two-time quantities such as the autocorrelation function C{t, s) or 
the response function R{t, s) only depend on the time difference t = t — s, where s (wait- 
ing time) is smaller than t (observation time), and both quantities are simply related to 
each other by the fluctuation-dissipation theorem 

fio,W = -i!5%M, (1.1) 

In the low-temperature phase (T < T^) the system undergoes phase ordering. In this non- 
equilibrium situation, C{t, s) and R{t, s) are non-trivial functions of both time variables, 
which only depend on their ratio at late times, i.e., in the self-similar domain growth (or 
coarsening) regime This behavior is usually referred to as aging 0. Moreover, no 
such simple relation as eq. (|1.1| ) holds between correlation and response, i.e., R(t, s) and 
dC{t, s)/ds are no longer proportional. It is then natural to characterize the distance to 
equilibrium of an aging system by the so-called fluctuation-dissipation ratio [0, ^ || 

TRit,s) 

dCjt, s) • ^ ' 

ds 

In recent years, several works p|, |, ||, |^, ||, 0, |, ||, |10| have been devoted to the study of the 
fluctuation-dissipation ratio for systems exhibiting domain growth, or for aging systems 
such as glasses and spin glasses, showing that in the low-temperature phase X(t, s) turns 
out to be a non-trivial function of its two arguments. In particular, for domain-growth 
systems, analytical and numerical studies indicate that the limit fluctuation-dissipation 
ratio, 

Xoo = lim lim X(t,s), (1.3) 

vanishes throughout the low-temperature phase ^, H- 

However, to date, only very little attention has been devoted to the response function 
R{t, s), and fluctuation-dissipation ratio X{t, s), for non-equilibrium systems at critical- 
ity. From now on, we will only have in mind ferromagnetic systems without quenched 
randomness. For instance one may wonder whether there exists, for a given model, a 
well-defined limit X^o at T = Tc, different from its trivial value X^o = in the low- 
temperature phase, and to what extent X^o is universal. Indeed a priori, for a system 
such as a ferromagnet, quenched from infinitely high temperature to its critical point. 
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the limit fluctuation-dissipation ratio X^o at T = (if it exists) may take any value 
between X^o = 1 (T > Tc. equilibrium) and Xoo = (T < Tc. domain growth). 

The only cases of critical systems for which the fluctuation-dissipation ratio has been 
considered are, to our knowledge, the models of ref. (random walk, free Gaussian field, 
and two-dimensional X-Y model at zero temperature) which share the limit fluctuation- 
dissipation ratio = 1/2) and the backgammon model, a mean-field model for which 
Tc = 0, where it has been shown that X^o = 1, up to a large logarithmic correction, for 
both energy fluctuations and density fluctuations [|TI1 , |12|. 

In a recent companion paper |]T3[, we have determined the fluctuation-dissipation 
ratio X{t, s) for the Glauber-Ising chain, another model for which Tc = 0. In particular, 
the limit fluctuation-dissipation ratio was found to be X^q = 1/2 (see also ref. [l^). In 
the present work we investigate the non-equilibrium correlation and response functions 
and the associated fluctuation-dissipation ratio in generic ferromagnetic models at their 
critical point. We first present (in section 2) an analytical study of the spherical model in 
arbitrary dimension. We then turn (in section 3) to a scaling analysis of the generic case, 
and to numerical simulations on the two-dimensional Ising model. One salient outcome 
of these joint works is the realization that the limit fluctuation-dissipation ratio X^o is a 
novel universal characteristic of critical dynamics, intrinsically related to non-equilibrium 
initial situations. 

The present paper is written in a self-contained fashion. For the spherical model, 
though our main intention lies in the study of non-equilibrium dynamics at the critical 
point, we shall present the three situations T > T^, T = T^, and T < in parallel. 
The latter case has already been the subject of a number of investigations, for both 
correlation and response 0, Results on the scaling behavior of the two-time auto- 
correlation function at Tc can be found in ref. |]15| . For the two-dimensional Ising model, 
several numerical works have already been devoted to its non-equilibrium dynamics in 
the low-temperature phase, concerning both correlations |TB]| and response We will 
therefore restrict our numerical study to the dynamics at the critical point. 



2 The spherical model 



2.1 Langevin dynamics 



The ferromagnetic spherical model was introduced by Berlin and Kac [|T7|, as an attempt 
to simplify the Ising model. It is solvable in any dimension, yet possesses non-trivial 



critical properties |117| , pTSf . Consider a lattice of points of arbitrary dimension D, chosen 
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to be hypercubic for simplicity, with unit lattice spacing. The spins S'x, situated at the 
lattice vertices x, are real variables subject to the constraint 



T.Sl = N, (2.1) 

X 

where is the number of spins in the system. The Hamiltonian of the model reads 

n = -J2 s.Sy, (2.2) 

(x,y) 

where the sum runs over pairs of neighboring sites. 

Throughout the following, we assume that the system is homogeneous, i.e., invariant 
under spatial translations. This holds for a finite sample with periodic boundary con- 
ditions, and (at least formally) for the infinite lattice. We also assume that the initial 
state of the system at t = is the infinite-temperature equilibrium state. This state is 
fully disordered, in the sense that spins are uncorrelated. The dynamics of the system is 
given by the stochastic differential Langevin equation 

^ = j:sy-ms.+vAt). (2.3) 

The first term, where y(x) denotes the 2D first neighbors of the site x, is equal to the 
gradient —dTL/dSy^, while X{t) is a Lagrange multiplier ensuring the constraint ( |2.1|) , 
which we choose to parameterize as 

X{t) = 2D + z{t), (2.4) 

and ?7x(t) is a Gaussian white noise with correlation 

{Ut)vy{t')) = 2T5^,y5{t-t'). (2.5) 



Equation (|2.3| ) can be solved in Fourier space. Defining the spatial Fourier transform 
by the formulas 

d^q 



r (q) = E /x /x = / ^ r (q) e^'^-, (2.6) 

where 

(2.7) 



f d^q 


r-^ dgi _ _ 




/ (27r)^ ~ , 


l^n 2tx ' ' 


J-n 2tx 



is the normalized integral over the first Brillouin zone, we obtain 

55^(q,t) 



^ -[a;(q) + ^(t)]5Mq,t)+r^Mq,t), (2. 
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where 

D 



o;(q) = 2^(l-cosgj q^, (2.9) 

a=l ^ 

and 

(r/^(q, t)r/^(q', t')) = 2T (27r)^ 5^(q + q') 6{t - t'). (2.10) 
The solution to eq. ( ^.8| ) reads 

5^(q,t) = e-(^)*-^W (^5^(q,t = 0) + J^^^'^^'^+^'^'^Y dt,^ , (2.11) 



with 



Z(t) = (\iti)dti. (2.12) 

JO 



2.2 Equal-time correlation function 

Our first goal is to compute the equal-time correlation function 

C^-y{t) = {S^{t)Sy{t)), (2.13) 

which is a function of the separation x — y, by translational invariance. We have in 
particular 

Co{t) = {S^{t)') = 1, (2.14) 
because of the spherical constraint ( |2.1| ), and 

Cx(t = 0) = 5x,o, (2.15) 



reflecting the absence of correlations in the initial state. In eq. ( |2.13D , the brackets denote 



the average over the ensemble of infinite-temperature initial configurations and over the 
thermal histories (realizations of the noise). 

In Fourier space the equal-time correlation function is defined by 

(^^(q,t)5^(q',t)) = (27r)^<5^(q + q')C^(q,t). (2.16) 



Using the expression (|2.11|) , averaging it over the white noise ?7^(q, t) with variance given 



by eq. ( |2.10| ), and imposing the condition 

C^(q,t = 0) = l (2.17) 



implied by eq. ( |2.15D , we obtain 



C^(q,t) = e-2-(q)*-2^W (^1 + 2T jJ^*e2.('i)*i+2Z(t,) ^^^^ _ ^2.18) 
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At this point, we are naturally led to introduce two functions, /(t) and g(T, t), which 
play a central role in the following developments. 

The function f{t) is explicitly given by 



^^'^ =1-^ = (e-^*/o(4t))'' (8v^t)-^/^ (2.19) 

where 

/o(^) = /^e"=°^9 ^ {27iz)-^/^e' (2.20) 

J 271 z—*oo 

is the modified Bessel function. 
The function 

g{T,t) = e^^^'^ (2.21) 
is related to f(t) by the constraint (|2.14| ), namely 

/(& ^'"i- = WJ) + '^i* - = ^- P-^^' 

which yields a linear Volterra integral equation for g{T,t) namely 

g{T, t) = fit) + 2T f fit - t,)giT, t,) dt,. (2.23) 
Jo 

This equation can be solved using temporal Laplace transforms, denoted by 

f\p)= fit)e-^'dt. (2.24) 
Jo 



We obtain 



with 

r d^n 1 

The dependence of g^iT,p) on temperature appears explicitly in eq. ( p.25| ). 

We now present an analysis of the long-time behavior of the function giT,t), consid- 
ering successively the paramagnetic phase (T > Tc), the ferromagnetic phase (T < Tc), 
and the critical point (T = Tc). To do so, we shall extensively utilize eq. (|2.25| ). We 
therefore investigate first the function f^ip), as given in eq. (|2.26| ). This function has no 
closed-form expression, except in one and two dimensions: 

D = l: f\p) - ^ 



pip + 8) 
7r|p + 8| Vb + 8| 



^ = 2 : /"(P) = K ( , (2.27) 
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where K is the complete elhptic integral. Together with the definition ( p.9|) , eq. ( p.26|) 
implies that f^{p) is analytic in the complex p-plane cut along the real interval [—8D, 0]. 
The behavior of f^{p) in the vicinity of the branch point at p = can be analyzed 
heuristically as follows. The asymptotic behavior of f{t) given in eq. ( p.l9| ) suggests that 
its Laplace transform has a universal singular part: 



/sg(p) ^ (87r)-^/2r(l - D/2K/2-\ (2.28) 
while there also exists a regular part of the form 

fl:,^{p) = A,- A2P + Asp"" + ■ ■ ■ , (2.29) 

where 

are non-universal (lattice-dependent) numbers, given in terms of integrals which are 
convergent for D — 2A; > 0. For instance Ai only exists for D > 2, and so on. Equa- 
tions ( |2.28| ) and ( p.29| ) jointly determine the small-p behavior of f^{p), as a function of 
the dimensionality D: 



D <2 
2<D <A 
D>A 



f\p) ^ (87r)-^/2p(^ _ D/2)p-^^-'"^\ 

f\p) ^A,- (87r)-^/^|r(l - D/2)b^/^-\ 

/L(p) ^ - A^p. (2.31) 



These expressions can be justified by more systematic studies (see e.g. ref. [0): f^ij)) 
possesses an asymptotic expansion involving only powers of the form and p^/^-i+n^ 
for n = 0, 1, . . . Whenever D = 2, 4, ... is an even integer, the two sequences of exponents 
merge, giving rise to logarithmic corrections, which shall be discarded throughout the 
following. 

In low enough dimension {D < 2), f^{p) diverges as p ^ 0. As a consequence, for 
any finite temperature, g^(T,p) has a pole at some positive value of p, denoted by l/r^ 



eq; 



away from the cut of f^{p). Hence 

9{T,t)^^ e*/^-, (2.32) 



and therefore, as further analyzed below, the system relaxes exponentially fast to equi- 
librium, with a finite relaxation time Teq. The latter diverges as the zero-temperature 
phase transition is approached, as 

-2/(2-D) 



(2(8vr)-^/2r(l - D/2)Ty ' ' . (2.33) 
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In high enough dimension {D > 2), f^{p = 0) = Ai is finite, so that the pole of 
g^{T,p) hits the cut of f^{p) at p = at a finite critical temperature 

n^^^(f4:^^y. (2.34) 



2Ai \J (27r)^ c<;(q 
As T — i> T+, the relaxation time Teq diverges according to 

'2(87r)-^/2|r(^_^/2)|T^ 



2/{D~2) 



2< D <A: Ten ^ I = — = 



c 

2 



2 A T 

^>4: ^^Tf^- (2.35) 

Note that these equations can be recast into the form Tcq ~ (T — Tc)"'^^'', where z/ is the 
critical exponent of the correlation length, equal to 1/{D ~2) for 2 < D < A and to 1/2 
for D > 4 |T^, while is the dynamic critical exponent, equal to 2 in the present cascQ 



We now discuss the asjTiiptotic behavior of the function g{T,t) according to temper- 
ature. Throughout the following, we will assume that D > 2, so that the model has a 
ferromagnetic transition at a finite Tc, given by eq. ( |2.34| ). 



In the paramagnetic phase (T > Tc), g{T,t) still grows exponentially, according to 
eq. (PI). 



In the ferromagnetic phase (T < T^), a careful analysis of eq. ( |2.25| ) yields 

. (2.36) 

eq eq 



where the spontaneous magnetization Meq is given by 

< = 1-^- (2.37) 

c 

At the critical point (T = T^), we obtain 

2<D<A: gin,t) ^ {D - 2){M)^''-' sin[(D - 2)7r/2] ^ f , 

c 

D>4: <,(T..t),-,^^. (2.38) 



summary of the values of static and dynamical exponents appearing in this work is given in 
Table 1. 
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Finally, eqs. (|2]2|), (|2]2|), O), and yield the following identities: 



oo 1 

fit)dt = — 

2lc 



/•oo 1 

1^ me-"-'-dt = — {T>%), 

We are now in a position to discuss the temporal behavior of the equal-time correlation 
function in the different phases. Its expression ( p.l8|) in Fourier space reads 

g-2a;(q)t 

using the definition ( p.21|) of g{T,t). We shall consider in particular the dynamical 
susceptibility 

Xit) = ^ E(^o(t)5.(t)) = = (2.41) 

for which eq. ( p.40| ) yields 



G^(q, t) = ——- (1 + 2T e'-^^^'^g{T, dhj , (2.40) 



The asymptotic expressions ( p.32| ), ( p.36| ), and ( p.38|) of g(T,t) lead to the following 
predictions. 

• In the paramagnetic phase (T > T^), the correlation function converges exponen- 
tially fast to its equilibrium value, which has the Ornstein-Zernike form 



where the equilibrium correlation length is given by 



Q = 2req. (2.44) 



The corresponding value of the equilibrium susceptibility is Xeq = ^eq = 2req. 
Eq. ( |2.43| ) implies an exponential and isotropic fall-off of correlations, of the form 
Cx.eq ~ e~l^l/^'"i, at large distances and for ^gq large, i.e., T close enough to T^. 



In the ferromagnetic phase (T < Tc), using the third of the identities (|2.39| ), we 
obtain a scaling form for the correlation function, namely 

C^(q, t) ^ (87rt)^/2 e'^'''', (2.45) 



9 



or equivalently, 

Cx(t)^M2qe-"'/(^*\ (2.46) 

in the regime where x is large (i.e., q is small) and t is large. Both the Gaussian 
profile of the correlation function, and its scaling law involving one single diverging 
length scale 

L{t) ~ t^/\ (2.47) 

reflect the diffusive nature of the coarsening process. The growing length L{t) 
can be interpreted as the characteristic size of an ordered domain. The dynamical 
susceptibility, 

Xit) ^ ^(8v^t)^/^ (2.48) 
grows as x(^) ~ L{t)^, or else as the volume explored by a diffusive process. 
At the critical point (T = Tc), the equilibrium correlation function reads 

cJq(q) ^ ^, (2.49) 

C (2 50) 

'-'^m- 47r^/2 |x|^-2' ^ ' 

These limiting expressions are reached according to scaling laws of the form 

C^{t) ~ Cx,eq*(xVt), (2.51) 



with 



2 < D < 4 : $(x) = 2x e'^^^^"") 2^/^-2 

Jo 



'0 



D>A: <l>(x) = l-e-2^ 

1 

The second expression of eq. (p.51|) has the general scaling form for the equal-time 
correlation function (see eq. (|3.9| )), with the known value of the static exponent 



of correlations rj = for the spherical model |T^, and with Zc = 2, already found 
above. 

The dynamical susceptibility grows linearly with time, as xit) ~ ^'(0)^5 i-^-; 

2<D<4: x{t)^^^t, 

D>A : x{t) ~ 2it. (2.53) 
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2.3 Two-time correlation function 



We now consider the two-time correlation function 



C^.y{t,S) = {S^{t)Sy{s)), 



(2.54) 



with < s (waiting time) < t (observation time). Its Fourier transform C^(q, t, s) is 
defined as in eq. (|2.16| ). Using eq. ( |2.11|) , we obtain 



C^(q,t,s) 



g{T,t)g{T,s) 



1 + 2T r e^'^^'i^''g{T,ti)dti 
Jo 



or else 



C^(q,t,s) = C^(q,s)e 



(2.55) 



(2.56) 



using the expression ( p.40| ) for C^(q, s). 

In the following, we shall be mostly interested in the two-time autocorrelation function 

C{t,s)^Co{t,s) = {S^{t)S^{s)) = |^^C^(q,t,s), (2.57) 



for which eq. ( p.55| ) yields 

Cit,s)= ^ 



g{T,t)g{T,s) 



f 



t + s 



t + s 



tAg{T,ti)dt, 



(2.58) 



The autocorrelation with the initial state assumes the simpler form 

f{t/2) 



C{t,s = 0) 



(2.59) 



The asymptotic expressions (|2.19| ), ( p.32|) , ( |2.36| ), and (p.38|) of the functions f(t) and 
g(T,t) lead to the following predictions. 



• In the paramagnetic phase (T > Tc), as s — *■ cxo with r = t — s fixed, the system 
converges to its equilibrium state, where the correlation function only depends on r: 

/CO 
/(ri/2)e-^^/(2r.,)^^^_ (2.60) 

This equilibrium correlation function decreases exponentially to zero as e~'^^^'^'^'"^^ 
when r ^ oo. The initial value Ccq(r = 0) = 1 is ensured by the second identity 
of eq. ( p39D . 
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In the ferromagnetic phase (T < Tc), two regimes need to be considered. In the 
first regime (s ^ oo and r fixed, i.e., 1 ~ r ^ s), using again the identities (|2.39| ), 
we obtain 

C{s + T,s)^ + (1 - M4)C,q,e(r), (2.61) 

where we have set 

fin/2) dn. (2.62) 

This function, which corresponds to the T ^ Tc hmit of eq. ( p.60|) , decreases only 
algebraically to zero when r — > cxd, as 

as implied by eq. (|2.19|) . The first identity of eq. ( p.39|) ensures that Ceq,c(T = 0) = 
1. 

In the second regime, where s and t are simultaneously large (i.e., 1 ^ s ~ r), with 
arbitrary ratio 

x=- = l + ->l, (2.64) 
s s 

the correlation function obeys a scaling law of the form 

When X ^ 1, this expression behaves as 

C{t,s) ^ AM^^x-^/\ (2.66) 

which can be recast into 



Cit, s) ~ Mi ( ) , (2.67) 



where L{t) is the length scale defined in eq. ( p.47| ), and A is the autocorrelation 
exponent (see section 3), which is equal to D/2 in the present case, in agreement 
with the result found in the n — oo limit of the 0(n) model (see ref. |jl|, p. 386, 
and references therein). 

Between these two regimes, the correlation function takes a plateau value 

T 



gEA = Jim lim C{s + t,s) = M^^ = 1 - t^, (2.68) 



known as the Edwards- Anderson order parameter (see e.g. ref. p][]). 
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Hereafter we shall refer to the first regime (1 ~ r ^ s) as the stationary regime, and 
to the second one (1 ^ s ~ r) as the scaling (or aging) regime. In the former, the 
system becomes stationary, though without reaching thermal equilibrium, because 
the system is coarsening. In the latter regime, as said above, the system is aging. 
It is possible to match these two kinds of behavior, corresponding respectively to 
eq. ( |2.61| ) and ( |2.65D , into a single expression: 

D/A 

(t + s)~ 



C{t = s + T, s)^{l- MI^)C,^,,{t) + [ ) , (2.69) 



which is the sum of a term corresponding to the stationary contribution, and a 
term corresponding to the aging one. Let us finally recall that, in the context of 
glassy dynamics, in a low-temperature phase, the first regime, where C(t, s) > qea, 
is usually referred to as the /? regime, while the second one, where C{t, s) < q^A, is 
referred to as the a regime [Q. 

At the critical point (T = Tc), the same two regimes are to be considered. However 
their physical interpretation is slightly different, since the order parameter M^q 
vanishes, and symmetry between the phases is restored. 

In the first regime (1 ~ r ^ s), the system again becomes stationary, the autocor- 
relation function behaving as the T — limit of eq. ( p.60| ), that is 



Cis + T,s) ^ Ceq,c(r), (2.70) 
which decreases algebraically to zero when t oo (cf. eq. (|2.6^)). In the second 



regime (1 <^ s ~ r), the correlation function obeys a scaling law of the form 

C{t,s)^T,s-^''/^-^^ F{x), (2.71) 
where the scaling function F{x) reads 

2<D<A: F{x) = 4(47r)-^/^ x'-^'/^x - 1)^-^/2, 

D>A: F{x) = ^^^^3^ {{x - 1)1-^/2 - (x + 1)^-^/^) . (2.72) 

In this regime the system is still aging, in the sense that C{t, s) bears a dependence 
in both time variables. However, the scaling of expression (p.71|) is different from 
that found in the low-temperature phase (see eq. (|2.65|) ), which depends on the ratio 
X = t/s only. The presence in eq. ( p. 711) of an additional s-dependence through 
the factor can be interpreted as coming from the anomalous dimension of 
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the field S'x at T^. In the critical region one has indeed Meq ~ (T - T^)^ ~ C'^/^ 
Replacing ^eq by s^^^" imphes the replacement of M^^ by 5-2/^/^^2= ^ g~{D-2+r))/zo_ 
With 7] = and = 2, the factor s~(-^/2-i) is thus recovered. Note that the static 
hyperscaling relation 2/5/z/ = D — 2 + rj holds for D < 4, while it is violated for 
D> 4 (see Table 1). 

Two limiting regimes are of interest. First, for x — 1, i.e., 1 ^ r ^ s, eq. ( p^.71D 



matches eq. ( |2.63|) . Second, for x ^ 1, i.e., 1 <^ s <^ t, one gets 

F{x) ^ Bx-^^/'% (2.73) 

where the autocorrelation exponent Ac (see section 3) is equal to 3D/2 — 2 if 2 < 
D < A, and to D above four dimensions, in agreement with the result found in 



ref. p 



We also quote for later reference the scaling law of the derivative 

dC{t, s) 



ds 



r,s-^/2Fi(x), (2.74) 



with 



D — 2 

Fi(x) = —F{x) -xF'{x), (2.75) 



2 < < 4 : F,ix) = (4.)-"/^ (^-^^"J-^|^'/^r'^ -^""''(^ - 

D>4: Fi{x) = {47Ty/^[{x-iy/^ + {x + iy/^). (2.76) 

2.4 Two-time response function 

Suppose now that the system is subjected to a small magnetic field H^^it), depending on 
the site x and on time t > in an arbitrary fashion. This amounts to adding to the 
ferromagnetic Hamiltonian ( |2.2[ ) a time- dependent perturbation of the form 

6n{t) = -J2H^{t)S^{t). (2.77) 

X 

The dynamics of the model is now given by the modified Langevin equation 

'^^^ T.Sy-mS. + H^{t)+r]^{t). (2.78) 



dt , . 
y(x) 
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Causality and invariance under spatial translations imply that we have, to first order 
in the magnetic field i^x(i), 

{S.{t))= fdsJ2R^^y{t,s)Hy{s) + --- (2.79) 

Jo y 

This formula defines the two-time response function Ry^_y{t,s) of the model. A more 
formal definition reads 



R^-y{t,s) 



6Hy{s) 



(2.80) 



{/fx(i) = 0} 

The solution to eq. ( |2.78| ) reads, in Fourier space. 



5^(q,t) = e-^^'i)*-^^ ( 5^(q,t = 0) + / e' 



^<^(q)tl+Z(ti) 



if^(q,ti) +/(q,ti)J dti 

(2.81) 

It can be checked that the Lagrange function X{t), and hence z{t) and Z{t), remain 
unchanged, to first order in the magnetic field. As a consequence, the two-time response 
function reads, in Fourier transform, 

5(5^(q,t)) 



i?^(q,t, s) 



6H^q,s) 



-<^(q){i-s) 



9iT,s) 
\i9iT,t) 



(2.82) 



{H^{t)=0} 

(cf. eq. ( p.56|) ). In the following, we shall be mostly interested in the diagonal component 
of the response function, corresponding to coinciding points: 



R{t,s) = Ro{t, s) 



SHJs) 



With the notations (|2l9| ), ( ^ ), eq. (^^821) yields 



(2.83) 



R{t,s) = f 



t - s 



9iT,s) 



2 J\9iT,t)- 

The response function at zero waiting time assumes the simpler form (cf. eq. ( p.59D ) 

f{t/2) 



(2.84) 



R{t,s = 0) = C{t,s = 0) 



9{T,t) 



(2.85) 



The asymptotic expressions (|2.19| ), ( |2.32|) , ( |2.36| ), and ( |2.38|) of F{t) and g(T,t) lead 
to the following predictions. 

• In the paramagnetic phase (T > T^), at equilibrium, the response function only 
depends on r, according to 

i?eq(r)=/(r/2)e--/(2-.). (2.86) 

Moreover, it is related to the equilibrium correlation function Coq(T) of eq. ( [^.60D 
by the fiuctuation-dissipation theorem ( p..l|) , as it should. 
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In the ferromagnetic phase (T < T^), the two regimes defined in the previous 
section for the case of the autocorrelation function are still to be considered. In the 
stationary regime (1 ~ r ^ s), the response function behaves as the T — > Tc limit 
of eq. ( p.86D , namely 

i?e,,c(r) = /(r/2) = (2.87) 
so that the fiuctuation-dissipation theorem is valid. 

On the contrary, in the scaling regime (1 <^ s ~ t), the response function has the 
form 

R{t, s) ^ {An{t - s)y/\tlsf^ = (47rs)-^/2 _ ^)-d/2^d/4^ (2.88) 



which, when compared to the corresponding expression (|2.65| ) for the autocorre- 



lation function, demonstrates the violation of the fluctuation-dissipation theorem 
(see section 2.5). 

At the critical point (T = T^), in the stationary regime (1 ~ r ^ s), the response 
function still behaves as in eq. ( p.87| ), so that the fluctuation-dissipation theorem 
still holds. In the scaling regime (1 -C s ~ t), the response function obeys a scaling 
law of the form 

i?(t,s) ^ 5-^/2^2(2;), (2.89) 
where the scaling function F2{x) reads 

D>4 : F2{x) = (47r)-^/'(x - 1^^^ (2.90) 



Again two limiting regimes are of interest. For a; — >■ 1, the scaling result ( ^.901 ) 



matches eq. ( p. 8?] ). For a; ^ 1, one flnds the same power-law fall-off for the 
functions F{x), Fi{x), and -^2(0;), that is 

F{x) ~ Fi{x) ~ ^2(0;) ~ (2.91) 

(see sections 2.5 and 3). 



2.5 Fluctuation-dissipation ratio 

As already mentioned in the introduction, the violation of the fluctuation-dissipation 
theorem ( |1.1|) out of thermal equilibrium can be characterized by the fluctuation-dissi- 
pation ratio X{t, s), deflned in eq. ( |1.2| ). In the case of the spherical model, the results 
derived so far yield at once the following predictions. 
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In the paramagnetic phase (T > Tc), the system converges to an equihbrium state, 
where the fluctuation-dissipation theorem holds. In other words, the fluctuation- 
dissipation ratio converges toward its equilibrium value 

Xeq = 1. (2.92) 

In the ferromagnetic phase (T < T^), the fluctuation-dissipation theorem ( |1 . 1| ) is 
only valid in the stationary regime (1 ~ r ^ s). On the contrary, in the scaling 
regime (1 ^ s ~ r), the results (|2.65| ) and (|2.88| ) imply that the fluctuation-dissi- 
pation ratio falls off as 

In particular, the limit fluctuation-dissipation ratio introduced in eq. (|1.3| ) reads 

= 0. (2.94) 

At the critical point (T = Tc), the scaling laws ( ^.74[ ) and ( |2.89| ) imply that the 
fluctuation-dissipation ratio X{t,s) becomes asymptotically a smooth function of 
the time ratio x = t/s: 

X(t,s) ^ X(x) = ^^, (2.95) 

i.e., explicitly, 

2 < D < 4 : Xix) = — ^ 



2 

1 + 



D-2\x+l. 

D>A: Xix) = ,\-^^n/2 - (2-96) 

1 + 



.X + 1. 

The scaling law ( |2.95| ) interpolates between the equilibrium behavior 

X{X) ^ Xeq = 1 (2.97) 

+1 

in the stationary regime of relatively short time differences, and a non-trivial limit 
value 

X{x) Xoo (2.98) 

at large time differences, given by 

D-2 

2<D<A: X^ = — 

D>4: Xoo = ^. (2.99) 
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Further comments on the scahng behavior of the fluctuation-dissipation ratio will be 
made in section 3.2. 



3 The generic situation 
3.1 Aging below 

Let us first briefly sketch the description of the dynamical behavior of a ferromagnetic 
system quenched from a disordered initial state to a temperature T < Tc ^ |, ^ H. 

In the scaling regime (1 ^ s ~ t), the autocorrelation C{t,s) is expected to be a 
function of the ratio L(t)/L{s) only, where the length scale L{t) ~ t^/^ is the characteristic 
size of an ordered domain, and z is the growth exponent, equal to 2 for non-conserved 
dynamics. More precisely, 

C{t,s)=M^^f{t/s), (3.1) 

where the scaling function / is temperature independent. Furthermore we have, for 
X = t/s > 1, i.e., K s < t, 

f{x)^Ax~^/\ (3.2) 
where A is the autocorrelation exponent ||16|. For the spherical model, eqs. ( p.65|) 



and (12:671 ) match eqs. (|J) and (jSj), with X = D/2 and z = 2. 



As a consequence, we have 



dC{t, s) M, 



fi{x), (3.3) 



ds s 
with /i(x) = — x/'(x), so that, when x ^ 1, 

/l(x)^Ax-"/^ (3.4) 

with Ai = AX/z. 

Although the situation of the response R{t, s) is less clear-cut, it is however reasonable 
to make the scaling assumption 

i?(t,s)^s-i-"/2(x), (3.5) 
where a > is an unknown exponent, and again with the behavior 

f2{x) ^ A2X-^/' (3.6) 

when X ^ 1. 
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The scaling law (|3.5|) holds for the spherical model, with a = D/2 — 1, as can be 
seen from eq. (|2.88|) . Furthermore, for non-conserved dynamics, at least in the case of 



a discrete broken symmetry, like e.g. in the Ising model, it has been argued |g, ^ 
that the integrated response pit,s) (to be defined in eq. ( p.l9| )), scales as pit,s) ~ 
L{s)^^ (p{L(t) /L{s)). This corresponds to eq. with a = 1/z = 1/2. 

The scaling laws (|3.3|), ( ^.51 ) imply 

X{t,s)^s--g{x), (3.7) 



with g{x) = (T/Mgq) /2(a;)//i(a;), in agreement with eq. ( |2.93|) for the spherical model, 
and especially 

Xoo = 0. (3.8) 



3.2 Aging at 

Let us now turn to the situation where a ferromagnetic system is quenched from a dis- 
ordered initial state to its critical point. 

In such a circumstance, spatial correlations develop in the system, just as in the 
critical state, but only over a length scale which grows like t^^^", where Zc is the dynamic 
critical exponent. For example the equal-time correlation function has the scaling form 

C.(t) = 1x1-2/^/^^0(1x1/^1/^^), (3.9) 

where f3 and u are the usual static critical exponents. The scaling function 0(x) goes to 
a constant for x ^ 0, while it falls off exponentially to zero for x — > oo, i.e., on scales 
smaller than t^^^" the system looks critical, while on larger scales it is disordered. This 
behavior is illustrated in the case of the spherical model by eq. ( p.51|) , corresponding to 
2P/u = D - 2 and = 2 in eq. (|3j). 

In the scaling region of the two-time plane, where both times s and t are large and 
comparable (1 ^ s ~ t), with arbitrary ratio x = t/s, the two-time autocorrelation 
function C(t, s) is expected to obey a scaling law of the form (see the discussion below 
eq. (|T|), and ref. pf) 

C{t,s)^s-^^/'''^F{x). (3.10) 



When both time scales are well separated (1 <^ s <^ t, i.e., a; ^ 1), the scaling function 
F{x) falls off as 

F{x) ^ B x-^^/'% (3.11) 



where Ac is the critical autocorrelation exponent |2^, related to the (magnetization) 
initial-slip critical exponent 0c ||15|] by Ac = -D — ZcQc- 
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We thus have 



ds 

with Fi{x) = —{2l3/uZc)F{x) — xF'{x), so that, when x ^ 1, 



with 



Fi{x) ^ Bix-^^/'% (3.13) 



B,= '-^^^B. (3.14) 

l/Zc 



For the spherical model, eqs. (|2.71|) , ( p.73|) , and (p.74|) respectively match eqs. (|3.10| ), 
( PTT| )> and (PTT^ ), with = 3D/2 - 2 ii D < A, and = D ii D > A (see Table 1). 

The similarity between the results ( p.74| ) and (|2.89| ), obtained in the case of the 
spherical model, strongly suggests that dC{t, s)/ds and R{t,s) behave similarly in the 
generic case, i.e., one is lead to hypothesize that a scaling law of the form ( p.l2| ), with 
the same power-law fall-off ( p.l3|) , holds for the response, that is 

R{t,s)^^s-'-'^/^^^F,ix), (3.15) 

c 

with, when x ^ 1, 

F2ix) ^ B2x'^^/^\ (3.16) 

The scaling laws ( |3.12| ) and ( |3.15| ) imply then that the fluctuation-dissipation ratio 
only depends on the time ratio x throughout the scaling region: 

X{t,s)^X{x) = ^^, (3.17) 

where the scaling function X{x) is universal. It appears indeed as a dimensionless com- 
bination of scaling functions. In turn, the hypothesis (|3.16| ) implies that the limit fluc- 
tuation-dissipation ratio reads 

= X{oo) = 1^. (3.18) 

This number thus appears as a dimensionless amplitude ratio, in the usual sense of 
critical phenomena. It is therefore a novel universal quantity of non-equilibrium critical 
dynamics, as already claimed in ref. ||T3 |. 

In the case of the spherical model, the analytical treatment of section 2 corroborates 
the above analysis, and yields the quantitative predictions ( p.96| ) and ( |2.99|) . 

In order to perform a numerical evaluation of X^, one needs to measure the response. 
A convenient way to do so is to measure instead the dimensionless integrated response 
function 

p{t,s) = T r R{t,u)du. (3.19) 

^0 
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By eq. (|2.79|) , this quantity is proportional to the thermoremanent magnetization Mtrm, 
i.e., the magnetization of the system at time t obtained after applying a small magnetic 
field h, uniform and constant, between t = and t = s: 

MTKM{t,s)^^p{t,s). (3.20) 

The thermoremanent magnetization is a natural quantity to measure experimentally in 
spin glasses 0, and it is also accessible to numerical simulations, for systems with and 
without quenched randomness (see section 3.3). 

The scaling law ( ^.151 ) for the response function implies 

pit,s)^s-^^/'''^Fsix), (3.21) 
with F2{x) = —{2[3/uZc)F^{x) — xF^{x)^ so that, when x 3> 1, 

F3(x)^53X-"^/^^ (3.22) 



with 



A clear representation of the evolution of X(t, s) in time is provided by the parametric 
plot of pit, s) against C(t, s), obtained by varying t at fixed s 0, 0, H- For well-separated 
times in the scaling regime (i.e., 1 <^ s <^ t), the common power-law behavior ( |3.11| ), 
( P.13| ), ( |3.16|) , and ( |3.22|) implies that the limit fluctuation-dissipation ratio has the al- 
ternative expression 

^oo = ^, (3.24) 

which is equivalent to eq. ( |3.18| ), due to eqs. ( p.l4| ) and ( p.23| ). In other words, the 
relationship (|1.2| ) also holds in integral form, that is 

p(t,s) ^XooC(t,s), (3.25) 

in the regime 1 ^ s ^ t. The limit fluctuation-dissipation ratio can thus be measured 
as the slope of the parametric plot in the scaling region, i.e., near the origin of the C — p 
plane. Eq. ( p.25| ) is expected to hold as long as C and p are much smaller than the 
crossover scale 

C*(s) = C(2s, s) ~ 3-''^''''% (3.26) 

corresponding to r = s. This quantity provides a measure of the size of the critical 
region, giving thus a quantitative definition of the critical analogue of M^^, involved in 
the discussion below eq. (|2.72| ). 
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3.3 The two-dimensional Ising model: numerical simulations 



In order to check the vahdity of the scahng analysis made in the previous section, beyond 
the case of the spherical model, we have performed numerical simulations on the ferro- 
magnetic Ising model on the square lattice, evolving under heat-bath (Glauber) dynamics 
at its critical temperature Tc = 2/ ln(l + a/2) ^ 2.2692, starting from a disordered initial 
state. The rules of the dynamics are as follows. Consider a finite system, consisting 
oi N = L'^ spins = ±1 situated at the vertices x of a square lattice, with periodic 
boundary conditions. The Ising Hamiltonian reads 

H = - E ^x^y (3-27) 
(x.y) 

where the sum runs over pairs of neighboring sites. Heat-bath dynamics consists in 
updating the spins crx(t) according to the stochastic rule 

^ I ^ (3.28) 

[_1 withpr„b.l^i^=!^M)ZZa. 

where the local field /?.x(^) acting on a^it) reads 

h.it) = E ayit), (3.29) 
y(x) 

with y(x) denoting the four neighbors of site x. 

Let us give a brief summary of known facts on the dynamics of the Ising model. For 
T < Tc, numerical studies have shown that the scaling forms ( |3.1| ) and ( ^.2] ) hold, with 
z = 2 (non-conserved dynamics) and A ~ 1.25 [1^. The integrated response function 
(in another form, known as the ZFC magnetization) has been measured in ref. 0]. At 
T = Tc, the dynamic critical exponent reads Zc ~ 2.17 and the autocorrelation 
exponent ^ 1.59 p|, |25[. 

Our aim is now to verify the hypotheses made in section 3.2, especially the scaling 
laws ( p. 151 ) and ( |3.21| ) for the response function, and to demonstrate the existence of a 



non-trivial limit Xoo- 

Computing C{t, s) with good statistics is rather easy, while the computation of p(t, s) 
requires more effort. We have followed the lines of the method introduced in ref. 0]. In 
order to isolate the diagonal component of the response function, a quenched, spatially 
random magnetic field, is applied to the system from t = to t = s. This magnetic field 
is of the form iJx = ^o^x, with a constant small amplitude ho, and a quenched random 
modulation, £x = ±1 with equal probability, independently at each site x. The heat-bath 
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dynamical rules are modified by adding up the magnetic field H-^ to the local field hy^{t) 
of eq. dpI). We have then 

(£xfTx(t)) = ho r Rit, U)du=^ p{t, S) = MTRM(t, s), (3.30) 
JO I 

where the bar means an average with respect to the distribution of the modulation of 
the magnetic field. 
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Figure 1: Log- log plot of the critical autocorrelation function C{t, s) of the two-dimen- 
sional Ising model, against time ratio x = t/s, for several values of the waiting time s. 
Data are multiplied by g^^^'^^", in order to demonstrate collapse into the scaling function 
F{x) of eq. (|3.10| ). Straight line: exponent —X^jzc ~ —0.73 of the fall-off at large x. 

We have first checked the validity of the scaling laws ( p. 101) , and especially (|3.21| ). 
Figures 1 and 2 respectively show log-log plots of the autocorrelation function C (t, s) and 
of the corresponding integrated response function p(t, s), against the time ratio x = t/s, 
for several values of the waiting time s. For each value of s, the simulations are run 
up to t/s = 10, and data are averaged over at least 500 independent samples of size 
300 X 300. For the response function, the amplitude of the quenched magnetic field 
reads ho = 0.05. Multiplying the data by s^/s/i^^c^ ^j^]^ 2(3 /uzc ~ 0.115, gives good data 
collapse, thus producing a plot of the scaling functions F{x) and -^3(0;). The data follow 
a power-law fall-off at large values of x, with a slope in good agreement with the value 
—Xc/zc ~ —0.73, shown on the plots as a straight line. 

We then turned to an investigation of the parametric plot of these data in the C — p 
plane. At the qualitative level, this plot, shown in Figure 3 for several values of the 
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waiting time s, confirms our expectations. The stationary regime (1 ~ r ^ s, i.e., 
roughly speaking, C > C*{s)), corresponds to the right part of the plot. The symbols 
show the data for small integer values of the time difference, r = t — s = 0, . . . , 8, 
illustrating the fast decay of correlation and integrated response in the stationary regime. 
The rightmost points, corresponding to r = 0, i.e., C = C{s,s) = 1, are compatible with 
the scaling law 1 — p{s, s) ~ C*{s) ~ g-'^P/'^^c^ ^jj-^g validity of the fluctuation-dissipation 
theorem is testified by the unit slope of this part of the plot, shown as a full straight line. 
The aging regime (1 ^ s ~ t, i.e., roughly speaking, C < C*{s)), corresponds to the 
left part of the plot. As expected, the data crossover toward a non-trivial slope, equal to 
the limit fluctuation-dissipation ratio Xqo. The dashed line shows the slope Xoo = 0.26, 
obtained by the analysis described below. 

-0.5 I 1 < < < 1 




In X 

Figure 2: Log-log plot of the critical integrated response function p{t, s) of the two- 
dimensional Ising model, against time ratio x = t/s, for several values of the waiting 
time s. Data are multiplied by s^/s/i^zc^ order to demonstrate collapse into the scaling 
function F^^x) of eq. ( |3.21|) . Straight line: exponent —Xc/zc ~ —0.73 of the fall-off at 
large x. 

In order to obtain a quantitative prediction of the limit fluctuation-dissipation ratio 
Xoo, we have followed two approaches. Figure 4 depicts the local slope of the plot of 
Figure 3, i.e., the ratio p/C, against C, in the aging regime. The data for the largest 
available waiting time s = 200 have been discarded from the analysis because they appear 
as too noisy on that scale. The data look pretty linear all over the range presented in the 
plot. This precocious scaling is due to the fact that the exponent 2[3 /vZc ~ 0.115 is small. 
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Hence the size of the critical region, given by the estimate ( |3.26|) , is very large, at least 
for waiting times s accessible to computer simulations. We have indeed, for example, 
C*(100) = C(200, 100) ~ 0.24. The straight lines show a constrained least-square fit of 
the three series of data, imposing a common intercept. The value of this intercept yields 
the prediction X^o ~ 0.262. 



0.8 - 
0.6 - 
0.4 - 
0.2 - 

1 --'^ ' ' ' ' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1.0 

C 

Figure 3: Parametric plot of the integrated response p{t, s) against the autocorrelation 
C{t,s), using the data of Figures 1 and 2. Symbols: data for integer time differences 
r = t — s = 0, . . . , 8. Full line: unit slope corresponding to the fluctuation-dissipation 
theorem in the stationary regime. Dashed line: limit slope X^o = 0.26 (see text and 
Figures 4 and 5). 

We have also followed an alternative approach, aiming at subtracting most of the 
deviations of the ratio p/C with respect to its limit X^o at C — > 0. This can be done 
by incorporating the known limit of the stationary regime, i.e., p ~ 1 as C ^ 1, into 
a quadratic phenomenological formula: p ~ XooC + (1 — Xao)C'^- This formula can be 
rewritten as ~ (p-C^)/(C(l-C)), suggesting to plot (p-C^)/(C(l-C)) against C, 
instead of the mere ratio p/C. This has been done in Figure 5. As expected, the vertical 
scale has been considerably enlarged. In return this procedure increases the statistical 
noise on the data points. The straight lines again show a constrained least-square fit, 
yielding Xoo ~ 0.260. 

We can conclude from this numerical analysis that we have 

Xoo = 0.26 ± 0.01 (3.31) 
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for the ferromagnetic Ising model in two dimensions. 




4 Discussion 

In the present work we deah with the dynamics of ferromagnetic spin systems quenched 
from infinite temperature to their critical state. This study, exemplified by the exact 
analysis of the spherical model in any dimension D > 2, and by numerical simulations on 
the two-dimensional Ising model, complements that of the Glauber-Ising chain, presented 
in a companion paper |13[. The main results obtained in this work can be summarized 
as follows. 

In such a non-equilibrium situation, these systems are aging in the sense that their 
correlation and response functions depend non-trivially on the waiting time s as well 
as on the observation time t, whenever these two times are simultaneously large. The 
corresponding scaling laws (see eqs. ( p.lO|) , (|3.12|) , (|3.15|) , and (|3.21| )), involve powers of s, 



related to the static anomalous dimension of the magnetization, and universal scaling 
functions of the ratio x = t/s. In the regime of large time separations, i.e., 1 ^ s ^ t (or 
X ^ 1), these scaling functions fall off algebraically with the common exponent Xc/zc- 

The fluctuation-dissipation ratio X{t,s), characterizing the violation of the fluctua- 
tion-dissipation theorem, has a universal scaling form X{x), and, for well-separated times 
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in the aging regime, it assumes a limit value X^o equal to a dimensionless amplitude ratio 
(see eqs. ( p.l8| ) and (|3.24| )). Therefore, as announced in ref. [|^, X^o is a novel universal 
characteristic of critical dynamics, which is intrinsically related to the non-equilibrium 
initial condition of a critical quench from a disordered state. 
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Figure 5: Parametric plot of the combination (p — C^)/(C(1 — C)) against C. Straight 
lines: constrained least-square fit yielding X^o ~ 0.260. 

The ferromagnetic models studied in the present work turn out to have values of X^o 
in the range 



0<X^< 



(4.1) 



We have indeed X^ = 1 - 2/D if 2 < D < 4, and X^o = 1/2 for D > 4, for the spherical 
model, and Xoo ~ 0.26 ± 0.01 for the two-dimensional Ising model. Let us mention that 
preliminary simulations on the three-dimensional Ising model yield Xoo ~ 0.40. The 
backgammon, for which X^o = 1 JTT], |12|, thus belongs to another class of models. 

The mean-field value 

= i (4.2) 



obtained for the spherical model in dimension D > A, also holds for a variety of mod- 



els which are not mean- field-like, including the Glauber-Ising chain and the two- 
dimensional X-Y model at zero temperature 

Let us finally discuss a few open questions. It would be interesting to know whether 
there is an analogue for the present case of the results found for models with discontin- 
uous spin-glass transitions, where the violation of the fluctuation- dissipation theorem is 
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related to the configurational entropy |26]. One would also like to know the status of the 
quantity Xqo for non-equilibrium systems with quenched disorder, or for systems defined 
by dynamical rules without detailed balance. 

In principle the limit fiuctuation-dissipation ratio could be calculated by field- 
theoretical renormalization-group methods, generalizing the computations done for uni- 



versal amplitude ratios in usual static critical phenomena as series in either e = A—D 



or in 1/n for the n-component Heisenberg model, the spherical model corresponding to 
the n ^ oo limit. The dimensionless time ratio x = t/s, appearing in the two-time 
autocorrelation and response functions and fiuctuation-dissipation ratio, is a temporal 
analogue of aspect ratios, which play an important role in static critical phenomena and 
finite-size scaling theory p8|. One may therefore wonder whether the latter, and espe- 



cially its latest developments involving conformal and modular invariance, could be used 
in order to put constraints on non-equilibrium critical dynamics. Generalized symmetry 
groups, such as those introduced in ref. ||29|, may also play a role in this issue. 
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Table and caption 



exponent 


spherical {2 < D < A) 


spherical {D > A) 


Ising {D = 2) 


V 








1/4 


P 


1/2 


1/2 


1/8 


v 


l/{D-2) 


1/2 


1 


z 


2 


2 


2 


X 


D/2 


D/2 


^ 1.25 




2 


2 


^ 2.17 


Ac 


3/5/2-2 


D 


^ 1.59 


©c 


1-D/A 





« 0.19 



Table 1: Static and dynamical exponents of the ferromagnetic spherical model and of 
the two-dimensional Ising model. First group: usual static critical exponents rj, /3, and v 
(equilibrium). Second group: zero-temperature dynamical exponents z and A (coarsening 
below Tc). Third group: dynamic critical exponents Zc, Ac, and ©c (non-equilibrium 
critical dynamics). 
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